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Parker Instability in a Self-Gravitating Magnetized Gas Disk: 

I. Linear Stability Analysis 

Sang Min Lee^’"^ and S. S. Hong^’^’® 

ABSTRACT 

To be a formation mechanism of such large-scale structures as giant molec¬ 
ular clouds (GMCs) and HI superclouds, the classical Parker instability driven 
by external gravity has to overcome three major obstacles: The convective mo¬ 
tion accompanying the instability generates thin sheets than large condensations. 
The degree of density enhancement achieved by the instability is too low to make 
dense interstellar clouds. The time and the length scales of the instability are sig- 
nihcantly longer and larger than the estimated formation time and the observed 
mean separation of the GMCs, respectively. This paper examines whether a re¬ 
placement of the driving agent from the external to the self gravity might remove 
these obstacles by activating the gravitational instability in the Galactic ISM disk. 
Self-gravitating, magnetized, gas disk bound by a hot halo medium is subject to a 
Parker-type instability driven by the self gravity, the usual Jeans instability and 
the convection. If the disk is thicker than a certain critical value, the Parker-type 
instability due to perturbations of undular mode assists the Jeans gravitational 
instability to form large-scale structures. Under the external gravity, the convec¬ 
tion triggered by interchange mode perturbations is to grow most rapidly; under 
the self gravity, however, the Jeans gravitational instability can have a growth 
rate higher than the limiting growth rate of the convective instability. The self 
gravity can thus suppress the convective motions, and a cooperative action of 
the Jeans and the Parker instabilities can remove all the obstacles confronting 
the classical version of the Parker instability. The mass and mean separation of 
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the structures resulting from the odd-parity undular mode solution are shown to 
agree better with the HI superclouds than with the GMCs. We briefly discuss 
how inclusions of the external gravity and cosmic rays would modify behaviors 
of the odd-parity undular mode solution. 

Subject headings: Instabilities — ISM: giant molecular clouds — ISM: magnetic 
helds 


1. INTRODUCTION 

In seminal studies of Parker in mid 1960s, he proved that magnetized gas disk under the 
influence of a constant gravitational acceleration becomes unstable to long wavelength per¬ 
turbations (Parker 1966, 1967). The instability named after Parker was once thought to be 
the formation mechanism of giant molecular cloud complexes (GMCs) in the Galaxy (Parker 
1966; Mouschovias 1974; Appenzeller 1974; Mouschovias et al. 1974; Blitz & Shu 1980; Shi- 
bata & Matsumoto 1991; Handa et al. 1992). This notion was based mostly on the results 
of stability analysis done in two dimensions, where only undular mode can be taken into 
acconnt. Many snbseqnent analyses done in three dimensions (eg., Zweibel & Kulsrud 1975; 
Hnghes & Proctor 1988), however, called our attention to the interchange mode, which may 
qnickly drive the system convectively unstable. Since the convection grows faster for shorter 
wavelength pertnrbations, it tends to shred the Galactic ISM disk into thin sheets of low 
density before the undular mode of long wavelength perturbations collects enough material 
to form well-dehned large-scale structures near the midplane (Asseo et al. 1978; Lachieze- 
Rey et al. 1980). To accept the classical Parker instability as the formation mechanism of 
such galactic structures as the GMCs or the HI superclouds, one needs means to control the 
disruptive tendency of the interchange mode. The observed density in GMCs is abont a few 
tens of hydrogen molecules per cubic centimeter. On the other hand, three-dimensional MHD 
simnlations of the Parker instability in nniform and non-nniform external gravities render to 
the system only a factor 2 enhancement in density over its initial eqnilibrium valne, which is 
far too low to generate interstellar clouds of even moderate density (Kim et al. 1998; Franco 
et al. 2002). This low degree of density enhancement is one of the most serious obstacles 
the Parker instability should overcome to be the GMC formation mechanism. 

Under canonical conditions of the Galactic ISM disk(cf., Spitzer 1978), the classical 
Parker instability driven by a nniform external gravity has the minimum growth time of 
~10® years and the corresponding length scale ~1.6kpc, which are clearly incompatible 
with the GMC observations. It was in 1990’s that non-nniform models of gravity were 
tried ont as a more realistic driving agent than the nniform one (Giz & Shu 1993; Kim & 
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Hong 1998). By introducing a non-uniform external gravity, whose acceleration increases 
linearly with the vertical distance, ; 2 , from the midplane up to ~ 500 pc and stays constant 
there on, one may reduce the time and length scales nominally to the degree compatible 
with the GMCs (Hong & Kim 1997; Kim & Hong 1998). However, the reduction of scales 
was not enough to give candidacy of the GMC formation mechanism back to the classical 
Parker instability, because the reduced scale values are relevant for activities in the regions 
of high Galactic altitudes, where vertical acceleration of the chosen external gravity becomes 
strongest but the GMGs do not generally reside. 

The undular mode of perturbations triggers the Parker instability preferably in those 
disks, where magnetic held strength decreases rapidly with vertical distance, because it is 
easy for such disk to develop strong magnetic buoyancy (Parker 1955; Hughes & Proctor 
1988). Under the uniform external gravity the held and density decrease with exponentially; 
while under the non-uniform external gravity they follow roughly a Gaussian function of ; 2 . 
Therefore, it is easier to trigger the Parker instability under the non-uniform gravity than 
the uniform one. This is how a replacement of the external driving agent from the uniform to 
the non-uniform gravity has brought to the system a signihcant reduction in the time scale. 
However, the reduction of time scale doesn’t necessarily mean that dense condensations 
form near the midplane within a reasonably short period of time. It only means a quick 
development of dihuse chaotic features high above the central plane. 

This line of reasoning led us to replace the driving agent of the Parker instability from 
the external gravity originated mostly from the Galactic stars to the internal one of the ISM 
itself. We expect the self gravity will ease the problem of low degree of density enhancement 
the classical Parker instability faced. The self gravity would prevent the Galactic ISM disk 
from developing convective activities, while it would also activate the Jeans gravitational 
instability preferably in the disk midplane. In this paper we will examine time, length and 
mass scales involved with the Parker and the Jeans instabilities in a magnetized gas disk 
under the self gravity. 

There have been many studies on the effect of magnetic helds upon the Jeans gravita¬ 
tional instability. Ghandrasekhar & Fermi (1953) examined development of the gravitational 
instability in a uniformly magnetized, plane-parallel gas layer. They showed that the exis¬ 
tence of uniform helds doesn’t change the critical wavelength for onset of the Jeans instability 
but only reduces its growth rate. When non-uniform held is considered, the critical wave¬ 
length does change and the gas layer becomes more unstable than the uniformly magnetized 
case (Stodolkeiwicz 1963; Nakano 1988). Nakamura et ah (1991) showed that self-gravitating, 
magnetized gas disk is subject to the Parker-type and the Jeans instabilities simultaneously. 
Since all of these investigations were done in the plane dehned by the vertical and the mag- 
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netic field directions, they couldn’t address the question whether the self gravity could put 
a curb upon flaring motions of the convection. Elmegreen (1982) is the hrst who seriously 
examined possible roles of the self gravity in modifying the classical picture of the Parker 
instability. His stability analysis was, in principle, of the three-dimensional nature. However, 
he integrated the basic MHD equations over the whole extent of the disk, which made the 
analysis two-dimensional in reality. 

In a more recent study, Nagai et al. (1998) made a fully three-dimensional stability 
analysis onto a uniformly magnetized, self-gravitating gas disk. However, the Parker insta¬ 
bility could not be triggered in their systems, because the magnetic helds were not vertically 
stratihed. Instead they examined how the disk thickness might influence the alignment of 
hlamentary structures with respect to the magnetic held direction. Chou et al. (2000) carried 
out an extensive analysis on non-uniformly magnetized, self-gravitating gas layers, for some 
cases including the external gravity as well. Their linear analysis was further extended to 
nonlinear, three-dimensional MHD simulations. In their systems, however, the interchange 
mode didn’t have a chance to drive convective motions, because only isothermal perturba¬ 
tions were applied to the background equilibrium which was in the same isothermal condition. 
And contrary to our objectives their main concern was on parsec-scale structures in the in¬ 
terstellar medium. We took adiabatic perturbations in one of our earlier three-dimensional 
stability analyses and saw a possibility of controlling the convection by the self gravity (Lee 
& Hong 1999; Lee et al. 2004). 

In principle the Parker instability in the Galactic ISM disk should be studied under 
simultaneous inhuence of the external and the self gravities. In this study we take the 
self gravity as the sole driving agent of the instability and examine what modifications it 
could bring to the classical picture of the Parker instability. This paper is organized as 
follows: In § 2 we will construct an initial equilibrium configuration for the self-gravitating, 
magnetized, isothermal gas disk, and then perform linear stability analysis onto it. In §3 
we will interpret the resulting dispersion relations in terms of the convective, the Parker and 
the Jeans instabilities. In §4 we will assess whether introduction of the self gravity as the 
driving agent might ease the major problems confronting the classical Parker instability. In 
the same section the resulting scales of the maximum growth rate will be compared with the 
observations of HI superclouds and GMCs. In the last section we will conclude the paper 
with some discussions. Appendix A details on coefficients of the linearized perturbation 
equations, and in Appendix B stable wave modes are isolated. Marginal stability analysis is 
carried out in Appendix G. 
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2. BASIC FORMULATIONS 
2.1. Initial Equilibrium Model 


We model the local part of the Galactic ISM disk as a self-gravitating, magnetized, 
isothermal gas layer, which is initially in magnetohydrostatic equilibrium and bound by a 
hot halo at z = iLzs,. Unperturbed magnetic held i?o = [0, Bo{z), 0] has azimuthal component 
only and its strength varies with vertical distance. 


In such disk the vertical structure is described by the equation of magnetohydrostatic 
equilibrium, 

d 
dz 


Po(z) 


Bojzy 

8tt 


= -po{z) 


d^oiz) 

dz 


( 1 ) 


where P, y, and p denote gas pressure, self-gravitational potential, and gas density, respec¬ 
tively. The potential is to be obtained from the Poisson equation. 


d^ 

dz^ 


%Ijo{.z) = AnGpo^z) , 


and the gas pressure is related to density as 


( 2 ) 


^ o (^) = cIp^{z) 


(3) 


with the isothermal sound speed, Cs, taken to be equal to the rms velocity of interstellar 
‘coudlets.’ In this study cosmic rays are ignored; instead, we will focus on the role of self 
gravity. The ratio, a = B^^z)/8 tiPo{z), of magnetic to gas pressure is kept constant of 2 : for 
the initial equilibrium. Equations (1) through (3) then give us the equilibrium stratihcations 
of density and magnetic held 


Po{,z) = po(0) sec\y{z/H) , 

(4) 

Bo{z) = Po(0) sech(; 2 /P), 

(5) 

and the resulting gravitational potential 


ijo^z) = 47rGpo(0)P^ ln[cosh(; 2 /P)], 

(6) 

where scale height H is dehned by 


u 1 1 + a 

~''^V27rGpM' 

(7) 
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2.2. Linearization 

The linearized MHD equations for a self-gravitating magnetized gas disk are given by 

^ + V ■ {pov) = 0 , (8) 

dv 1 — 1 — — 

Po-^ = -VPi — (V X Pi) X Po + ^ (V X Po) X Pi - piV^/’o - PoV^/'i, (9) 

at All An 

no 

^ = Vx(TxPo), (10) 

BP 

+-fPoW ■ V + {v ■W)Po = Q , (11) 


V^^I = dTTGpi, (12) 

where all the symbols have their usual meanings and the ones with subscript 1 denote their 
perturbed quantities. We take adiabatic relation P = np^ as the equation of state for the 
perturbed material [Eq. (11)]. 

Seeking solutions of the form, qi{x, y, z, t) = qi{z) exp [i{kxX + kyp — ut)], we recast the 
linearized equations in dimensionless forms as 

hi—= + pUy + 2iQuz — i ) (13) 

po dC 

^ Pi 2ap Pi^ Pi^ 

flMx = (14) 

1 H" Q; Pq 1 -|- o; Pq 1 H- o; Pq 

p Pi .2aQBi, 

_ 1 d Pi 20 Pi . 2 ap Biz 2q; d Biy 

X \ I 'll =■ _ _ _ — _ _ — X _ _ -I- __— 

1 -|- o; Pq 1 -|- o; Pq 1 -[- q; Pq 1 -h o; (P^ Pq 

1 + a Bo dQ Po 

kl^ = -puz;, (16) 

Pq 

^ + i0M2 + , (17) 


Q^lZ 

= -PUZ: 


(18) 
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VL^ = 2i + l{iu^ + 'nuy) - *7^ , (19) 

and finally 

(y + 7)v>= 2(1-0=)^. (20) 

«C po 

Here we have taken 77, yjl + a Cg (= H) and VjH as length, velocity and frequency units, 
respectively: ^ = KH, 7] = kyH,C = z/H,Ux = v^/V, Uy = Vy/V, = v^/V, (p = 
Tpi/V"^, n = u{V/H)~^, and 0 = tanh(^. Eliminating all the variables in favor of Uz and 
p, we combine equations (13) through (20) into the following two second-order ordinary 
differential equations: 


d 


AqUz = Bi-^p + BoP 

dC 


( 21 ) 


and ^ 

+ CoP = Di—Uz + DqUz- ( 22 ) 

ctC dC 

Coefficients Aq, Ai, A 2 , Bq, Bi, Cq, C 2 , Dq, and Di are given in Appendix A. 


2.3. Boundary Conditions and Numerical Method 


To £x boundary conditions we will follow Goldreich & Linden-Bell (1965) and Tomisaka 
& Ikeuchi (1983). Let us suppose that the upper boundary of the disk is displaced from 2 : = 
±Za by a small amount 6z = s exp [i {k^x + kytj — ut)] with s H, which will then give us 
the vertical component of perturbation velocity at the boundary as 

Vz{Za) = -iuj Sz. (23) 


An application of the Gauss flux theorem to the perturbed boundary gives us 

d'lfji 


-kl/jiiZa) - 


dz 


= 4:7rGpo{za)6z, 


(24) 


since 'ipi{za+) satisfies the Laplace equation in the rarefied halo. 


A third condition comes from pressure continuity. It is the sum of the gas and the 
magnetic pressures that should be kept continuous at the boundary. By integrating 2 ;- 
component of the MHD momentum equation from Za to Za + Sz, one may have the total 
pressure at the perturbed boundary as 


P{za + Sz) + 

on 


\za+Sz 


P(z.) + 


Po{,Za) 


dipo 

dz 


Sz + 0{Sz^). 



Therefore, the total boundary pressure, Ptot^, calculated from the disk side becomes 


P,TiZa + fc) = PoiZa) + PliZa) + —B, 

8n 


+ - BoBly 

dTT ^ 


- Po{Za) 


dlpo 


dz 


Sz 


( 25 ) 


to hrst order of 5z. For the same boundary pressure, °, but calculated from the halo 
side, we may have 


P 


+ dz) — Po{Za) + 


halo 




(26) 


where perturbation of the halo gas pressure Pi^^°{za) has been ignored, as the halo medium 
is thought to be very hot. Also ignored is the gravitational potential term, because the halo 
is hlled with a rarehed medium. By equating P^^^ to P^° at the boundary, we hnally have 
the pressure continuity condition as 


PliZa 


An 


-B^B 


Q-D\y 


- Po{Za) 


d'lpo 


dz 




(27) 


Now B\y'^'^ has to be specihed. From the induction equation [Eq. (10)] we obtain 


B\ VxBq Cx 

LO 


— ^ {v^Bo) - — v^Bo 

OJ OZ 00 


00 


(Oy VzBf) C2 


(28) 


On the other hand, no current density is expected in the halo, because matter density there 
is extremely low. Therefore, can be given by magnetic potential 0^ simply as 

(29) 


with 


0m = A exp {ikxX + ikyp — kz) and k'^ = kl + ky. 


(30) 


Equation (28) tells us that the vertical component of the magnetic field becomes Bff'^ = 
— {ky/00) VzBo at the disk side boundary. From the halo side. Equation (29) shows B\f° = 
—kBo(f)m- By equating Bff^ to B'^f°, we have 0m = {ky/koo)vz and the azimuthal held 
component B\^‘^ = i{ky/ku) VzBq. And its substitution for B^y^° in equation (27) hnally 
gives the pressure continuity condition as: 


PliZa 


— BoBly 

An 


Po{za 


#0 


dz 


5z = 

X An 


A 

koo 


VAZa 


(31) 


In terms of the dimensionless variables the upper boundary conditions [Eqs. (23), (24), 
and (31)] take the following forms: 


iaSC, = !l,(Ca), 


( 32 ) 
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and 


-7 iil V cosh^C( 


-k^fiCa 


(Pif 


dip 

dC 


= 2sech^{Ca)SC, 


dC 


iVt - 7 "Drcosh Ca - 2ai^Vt^ (p((a) 


- 2 (1 + a)einvsc + 2a{n^- 


2a 


1 + a J d( 


20(7 -1)V + - 20a <! - 

fz 


7 - 2^2 , 2 a+ 7 


1 -|- Q; 


r + 


1 “h q; 


du 

Yc 

2 a 7 


Ca 


2 1 n 2 I I 2 ; 2 

" l« +7T^'''= 


where 


zjH = Ca and V = 


2q: + 7 
1 “h O; 


k^n^ + 


2 a 7 

(1 + a )2 




(33) 


Uz{Ca) 0 , 

(34) 


At the midplane, parity conditions shonld be specihed: If symmetric (odd-parity) soln 
tion is songht for vertical velocity Uz, it shonld satisfy 

dip 


dC 


C=o 


= 0 and Uz{C = 0 ) = 0 


at the midplane; while anti-symmetric (even-parity) one reqnires 

duz 


ip{C = 0 ) = 0 and 


dC 


= 0 . 


C=o 


(35) 


(36) 


The odd-parity solntions are mirror symmetric with respect to the midplane, and the even- 
parity ones are of the midplane crossing type. 


Along with the npper-bonndary continnity and the midplane parity conditions, eqna- 
tions (21) and (22) comprise an eigenvalne problem for For given perturbation wave 
vector (C, 7 ) we take a trial value for and integrate the equations from C = Ca to C = 
0 with the fourth-order Runge-Kutta method. If the resulting solution satishes the chosen 
parity condition at 2 : = 0, the trial is adopted as the eigenvalue. If not, we take another trial 
and perform the same integration until the solution fulhlls the midplane parity condition. 

The resulting growth rate squared, —is illustrated as a function of Ux = CfanhCa 
and Uy = T^tanhCa (not of C and 7 ), which will keep the slope of dispersion relation the 
same in limits Ux ^ 0 or Uy ^ 0, regardless of the disk thickness (Elmegreen & Elmegreen 
1978). If we notice total column density is 2po(0) x H for unbound isothermal disk and 
2po(0) X H tanhCa for halo-bound one, we may call Ux and Uy normalized wave numbers that 
are effective for the halo-bound disk. The three parameters (a, 7 , Ca) characterize the disk 
system and the normalized effective wave numbers (z/j,, Uy) do the perturbation. 
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3. CONVECTION, PARKER AND JEANS INSTABILITIES 

Whenever a need arises to differentiate the undular-mode instability driven by external 
gravity from that by self gravity, we will use the term classical Parker, for which both odd- 
and even-parities are applicable. In the context of self gravity, we will call the instability 
triggered by the same undular mode but even-parity perturbations as simply the Parker or 
the pure Parker instability, while reserving the term Parker-Jeans for the odd-parity undular 
mode. On the other hand, the instability triggered by short wavelength perturbations of 
the interchange mode will be called broadly as convection. Long wavelength perturbations 
of the same mode may of course activate the gravitational instability. In this section we 
will examine how the convection, the pure Parker instability, and the Jeans gravitational 
instability compete or collaborate with each other in forming or de-structuring large-scale 
structures. 


3.1. Undular Mode 

To have a comparison standard we first calculate the dispersion relation for the pure 
gravitational instability by taking a = 0, = 0 and imposing the odd-parity condition. For 

several cases of disk thickness. Figure la illustrates how the normalized growth rate of the 
gravitational instability changes with the azimuthal effective wave number. In general, the 
Jeans instability grows faster in thicker disk. However, once it is thicker than about 4 times 
the scale height, the growth rate becomes almost independent of the disk thickness. 

With the same set of system parameters except for a = 1.0, we calculate —again 
as a function of Uy. The resulting dispersion relations for odd- and even-parity solutions 
are shown, in Figure lb, by solid and dotted lines, respectively. Because even-parity cannot 
activate the gravitational instability (Simon 1965; Elmegreen & Elmegreen 1978), by simply 
changing the parity condition from even to odd, we may easily differentiate role of the Parker 
instability from that of the Jeans in the dispersion relation. The dotted lines in Figure lb 
are, therefore, for the Parker instability only; while the solid ones are for the case where 
both the Parker and the Jeans instabilities operate simultaneously. This led us to coin the 
term Parker-Jeans especially for the odd-parity undular mode. 

A comparison of the solid lines in Figures la and lb indicates that the presence of 
magnetic helds in moderately thick disks boosts the growth rate of the gravitational insta¬ 
bility signihcantly. The Parker and the Jeans instabilities interact constructively in the disks 
thicker than, say, one half the scale height; while in the disks thinner than that, the same 
magnetic helds now add rigidity or extra pressure to the gas, thereby hindering growth of 
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the gravitational instability. With marginal stability analysis (see Appendix C) we have 
determined the critical thickness, over which the two instabilities interact constructively or 
destructively: It is 0.549 for a moderately magnetized {a ~ 1), isothermal (7 = 1) disk. 
As can be seen from the dotted lines in Figure lb, growth rate of the Parker instability 
depends on thickness more sensitively than that of the Jeans instability. The dotted line 
almost coinciding with the horizontal axis is for = 1 , and the dispersion relation of even- 
parity becomes hardly noticeable for (a < 1-0. Comparison of the dotted (even-parity) and 
solid (odd-parity) lines in Figure lb suggests that a marked increase can be achieved in the 
growth rate of the Parker instability by triggering the gravitational instability. This opens 
up a possibility of alleviating the time scale problem of the classical Parker instability. 

To understand how the magnetic held strength modihes the dispersion relation for the 
undular mode, we hx the disk thickness at 5iJ and vary a over 0, 1, 2, 3, and 4. The 
resulting —is plotted, in Figure 2 , against the effective azimuthal wave number. As 
before the solid and dotted lines are for odd- and even-parity solutions, respectively. It 
is clear from the hgure that the maximum growth rate of the Parker instability increases 
with held strength (Notice diherence in the ordinates between Figures 1 and 2.). Under 
an inhuence of vertical acceleration, the undulated held lines give the magnetic buoyancy a 
chance to unload the gas material oh their shoulders towards the midplane valleys. As the 
held strength increases, so does the buoyancy, and hence the system becomes highly unstable 
to the perturbations of undular mode. The resulting escalation of the Parker instability then 
boosts development of the gravitational instability. The maximum growth rate of the odd- 
parity solutions, therefore, increases signihcantly with a. The critical wave number for 
marginal stability is known to increase with a as (1 -|- 2aY^‘^ (Stodolkiewicz 1963), which 
also indicates the constructive nature of Parker and Jeans interaction. 

We now turn our attention to the remaining system parameter, 7 , which was intro¬ 
duced as a measure of material rigidity. To trigger the classical Parker instability under 
external gravity, effective ratio of the specihc heats for the perturbed material should be 
less than a certain critical value, 7 c. Many authors have shown 7 c = 1 -|- a under external 
gravity (Newcomb 1961; Parker 1966; Gilman 1970; Kim & Hong 1998). In a recent study 
Lee (2002) proved that the critical value takes the same 1 + a under the influence of self 
gravity as well. To somewhat dramatize the effect of 7 on the growth rate we decided to 
hx q; at a rather small value 0.1 and examined the dispersion relation with 7 varying 0.68 
through 1.1. The solid and dotted lines in Figure 3 are for the odd- and even-parity solutions, 
respectively. For soft materials the dispersion curve exhibits two maxima: The gravitational 
instability is responsible for the peak occurring at smaller wave number and the Parker is 
for the other at larger one. As 7 increases, height of the Parker peak decreases. This is 
because, as the material becomes harder, the magnetized gas system becomes more stable to 
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the undular mode perturbation. If 7 exceeds its critical value, the Parker instability cannot 
be triggered at all. Owing to a heightening activity of the Parker instability in softer media, 
the Jeans peak becomes taller as 7 decreases. 

For most cases shown in Figure 3, the dotted lines of even-parity match smoothly with 
the solid lines of odd-parity. For somewhat large 7 ’s, however, the even-parity perturbations 
of undular mode grow a little faster than the odd-parity ones. Under self gravity the growth 
rate difference between the even- and odd-parity solutions is too small to render any signif¬ 
icant consequences to the spatial distribution of the resulting condensations with respect to 
the midplane. In the case of classical Parker instability, however, the even-parity perturba¬ 
tions grow substantially faster than the odd-parity ones. Therefore, the condensations, if 
they were generated by the classical Parker instability, would place their maximum density 
points on the northern and southern sides of the midplane alternatively (Kim et ah 2000). 

Figure 3 suggests that the undular mode of odd-parity triggers both the Jeans and 
the Parker instabilities in the 7 = 0.68 disk, for example. The two most rapidly growing 
perturbations in the disk have ~ 8.6 H for the Jeans instability and Xy ~ 3.5 H for the 
Parker case. On the other hand, the perturbations of the same undular mode with the 
same Xy ~ 8.6 H but of even-parity would turn off the gravitational instability and drive the 
system to develop only the Parker. We have shown, in Figures 4a, 4b and 4c, eigen-solutions 
for the cases with Xy ~ 8.6 H of odd-parity, with Xy ~ 3.5 H of odd-parity, and with Xy ~ 
8.6 H of even-parity, respectively. In the left panel, the color changes from dark red to white 
as density increases, and the numbers marked on each contours represent logarithm of the 
density perturbation, pi{y, z), measured in units of the initial equilibrium density, po{z = 0 ), 
at the midplane. In the right hand panels, contours in thin solid lines are for the magnetic 
fields and the ones in thick lines delineate the halo boundary. Small arrows inside the panels 
represent velocity vectors, a big arrow on top of frame a and b corresponds to two times the 
sound speed, and that of frame c does just one. 

An action of the self gravity is apparent from the velocity pattern in Figure 4a. In 
the region of inflated held lines the magnetic buoyancy drives matter to rise against the 
self gravity; while in the magnetic valley the self gravity does the driving in an opposite 
direction and collects material towards the center. We see motions taking place everywhere 
in the yz-plane. Even in the midplane one may witness a converging motion, which is made 
possible by the self gravity and cannot be expected under the external gravity. As a hnal 
product of the Jeans instability the system will bear out a dense clump centered at Y ~ 4.3. 
Velocity vectors in Figure 4b exhibit a pattern markedly different from the ones in Figure 4a. 
Magnitude of the velocity vectors quickly diminishes towards central plane, and they hardly 
exhibit any motions in the midplane. The magnetic held lines in low altitudes maintain 
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their initial configuration almost intact. Features of the perturbed density, the magnetic 
fields and the velocity patterns in Figure 4b are all very similar to the corresponding cases 
of the classical Parker instability that is driven by an external gravity of linear model (Lee 
1993). As a final product of the Parker instability the system will bear out clumps centered 
at Y ~ 1.7, 5.2 and 8.5. The iso-density contours for odd-parity solutions in Figures 4a and 
4b demonstrate dense clumps formed right in the midplane. The deformation of boundary is 
closely related to the enhanced density near the midplane (Elmegreen & Elmegreen 1978). A 
rise and a fall of the boundary surface from its initial equilibrium height Za = 5 H accompany 
a fall and a rise in the midplane density, respectively. 

Eigen-solutions of the even-parity shown in Figure 4c make a marked contrast to the 
odd-parity solutions in Figures 4a and 4b. The upper and lower boundaries make the rise- 
and-fall motions in unison. Active motions are apparent only in high above the midplane. 
These features are very much reminiscent of the Parker instability driven by a linear model 
of external gravity (Lee 1993). It is quite clear that the classical Parker instability may not 
be the formation mechanism of such large-scale structures as GMCs or HI superclouds. The 
GMCs are attributes of the Galactic midplane, while the Parker instability generates actions 
at high altitudes. To have such dense structures as GMGs in the disk central plane, we need 
to invoke the self gravity not only to drive the Parker undular mode but also to trigger the 
Jeans gravitational instability. 

As far as forming large-scale structures is concerned, the odd-parity undular mode de¬ 
serves our special attention, because it will drive the system to develop the Jeans gravitational 
instability and the Parker instability simultaneously. If the two instabilities work construc¬ 
tively, they would generate dense condensations in the central plane. On the other hand, 
the undular mode perturbations of even-parity drive the system to develop only the Parker 
instability, whose activities are confined mostly to high altitude regions, where GMGs ar not 
generally found. The question then is whether the resulting condensations of the odd-parity 
solutions are safe from disruption by convective motions. 


3.2. Interchange Mode 

To construct dispersion relation for the pure interchange mode, we fix azimuthal wave 
number at zero {vy = 0) and calculate the growth rate as a function of effective radial wave 
number The results for the odd-parity solutions are shown, in Figure 5, for the disks of a 
moderately magnetized (a = 1), self-gravitating, isothermal gas medium. As was seen from 
Figure la for the non-magnetized gas disk, the maximum growth rate generally increases 
with disk thickness. However, for the disks thicker than ~5 H, the dispersion relations for 



different thickness are hardly distinguishable from each other. Once the upper boundary 
is placed beyond a certain distance from the midplane, development of the gravitational 
instability in the densest part of the system doesn’t depend much on where the boundary 
is located. In disks thinner than the growth rate again becomes independent of 

the thickness. In such thin halo-bound disks, there is practically no gradient in density 
over vertical distance, and the gas acts like an incompressible fluid (Goldreich & Lynden- 
Bell 1965; Elmegreen & Elmegreen 1978; Lubow & Pringle 1993; Nagai et ah 1998; Lee 
2002). A sizeable difference exists in the maximum growth rate between non-magnetized 
(see Figure la) and magnetized (see Figure 5) disks. Perturbations of the interchange mode 
make density and magnetic flux change in phase, and thus activate fast magneto-acoustic 
waves to propagate in the radial direction. Since the fast waves tend to disperse enhanced 
density, the gravitational instability triggered by the interchange mode ought to grow slower 
than that by the undular mode. This is the reason for the difference between the two modes. 

We saw from Figure 2 that the existence of magnetic fields boosts the gravitational 
instability triggered by the undular mode. To check whether the same boosting is possible 
for the interchange mode, we calculate —of the odd-parity solutions for a = 0 , 1 , 2 , and 
3. As shown in Figure 6 , the maximum growth rate decreases with a. This trend is opposite 
to the case of the undular mode. With the interchange perturbations the magnetic helds 
have no means for adding extra material to the densest region of the disk; instead they act 
to disperse the enhanced density by driving the fast magneto-acoustic waves. Therefore, the 
growth rate of the interchange mode decreases with increasing held strength. 

In all the cases illustrated in Figures 5 and 6 , the magnetic held is too strong to satisfy 
the criterion, 7 < 1 — o;, for convection. Originally this criterion is known to hold for disks 
under external gravity (Newcomb 1961; Parker 1967; Zweibel & Kulsrud 1975; Asseo et al. 
1978; Lachieze-Rey et al. 1978). As for the 7 c criterion for the classical Parker instability, 
Lee (2002) showed that the same convection criterion, 7 < 1 — a, is applicable to the self 
gravity case. 

Adopting a weakly magnetized {a = 0.1) disk to see the convective motions and assigning 
small values to the ratio of ehective specihc heats to meet the convection criterion, we 
calculated the dispersion relation as a function of The results are shown in Figure 7, 
where the solid and the dotted lines are for the odd- and the even-parities, respectively, 
and on each curve 7 is marked. The well-dehned peak corresponds to the Jeans instability; 
while the monotonically increasing part is owing to the convective instability. The even- 
parity perturbations cannot trigger the gravitational instability but certainly can activate 
the convection. After passing a minimum the growth rate hrst increases monotonically with 
radial wave number and reaches a hnite value asymptotically as the wave number further 



increases: Arrows outside the frame indicate the limiting values of the growth rate. Under 
external gravities the interchange mode perturbations of the shortest wavelengths are to 
grow most rapidly and deprive the undular mode of the chance to form large-scale structures. 
Under the self gravity, however, the maximum growth rate of the gravitational instability 
can be higher than the limiting growth rate of the convective instability. The self gravity 
may thus have a chance fully to suppress convective motions. This is a most important 
attribute of the self-gravitating disk. Since the external gravity exerts acceleration of hxed 
strength only in the vertical direction, there is no means of de-structuring the vertical sheets 
generated by the convective motions. Once a substantial degree of density enhancement is 
achieved, however, the ever increasing self gravity can provide the system with a means for 
redirecting the vertical motions towards the dense mass center. 

Even if the convection can grow faster than the Jeans instability in linear stage of 
the disk evolution, particularly in ‘soft’ medium of small 7 ’s, the convective motion won’t 
prevent the system from forming dense structures in the midplane for two reasons: First, 
no buoyancy force is available in the midplane and hence no convective motions can take 
place there. Second, the acceleration due to self gravity increases with density enhancement. 
Therefore, once the gravitational instability is triggered by odd-parity perturbations, it will 
eventually dominate the disk over the convection. 

The eigen-solution of the odd-parity interchange mode was constructed for the case with 
7 = 0.68 and = 0.6 (A^, = 10.5 iJ), which corresponds the well-dehned growth rate peak 
in Figure 7. The same was done of the even-parity interchange mode with = 3.0 (A^, 
= 2.177), which is a convection dominating case. Figures 8 a and 8 b compare the two sets 
of eigen-solutions in the xz-plane. The left and right panels illustrate density and velocity 
fields, respectively, and the color changes from dark red to white as density increases. The 
numbers marked to thin solid lines on the left are logarithms of pi{x,z)/po{z = 0). The 
thick solid lines on the right demarcate the disk-halo boundary, the arrows inside the frames 
represent velocity vectors, and the horizontal arrow on top of each frame is a measure of the 
sound speed. 

A well-dehned density clump in Figure 8 a is obviously a result of the converging motions 
driven by the Jeans instability. The velocity pattern in Figure 8 a closely resembles the same 
pattern in Figure 4a: Converging motions are evident all over the xz-plaxie in Figure 8 a and 
the yz-p\a.ne in Figure 4a. Figure 8 a suggests that once the interchange mode perturbation 
generates a slight density enhancement at a midplane point, towards which the self gravity 
would collect material and form a density clump. On the other hand, we don’t expect from 
Figure 8 b such clump to form in the midplane: Motions are absent from the low altitude 
region. 



3.3. Mixed Mode 


We now calculate the growth rate as a function of Vx and Vy. The results are presented 
in a three dimensional surface plot of —over the wave number domain. As our interest is 
primarily in forming large-scale structures in the Galactic ISM disk, we £x the disk thickness 
at Ca = 5, which is sufficiently larger than the critical value, 0.549 (see §3.1). Figures 9a and 
9b are for even-parity solutions with (a = 1, 7 = 1) and (a = 0.1, 7 = 0.8), respectively. 
The surface plot in Figure 9a looks like a long ‘tunnel’, and we already met, in Figure 
lb, the dispersion curve outlining the tunnel ‘entrance’ at Vx = 0. Cross-section of the 
tunnel cut perpendicularly to the z/a;-axis describes the dispersion relation of the Parker 
instability triggered by undular portion of the mixed mode perturbation. Please note that 
the tunnel ridge inclines slowly with Vx- This is owing to the interchange part of the mixed 
mode. The pure interchange mode (i.e., Vy = 0) doesn’t trigger the magnetic Rayleigh- 
Taylor instability (cf., Chandrasekhar [1961] for this terminology), since this disk does not 
satisfy the 7 < 1 — a criterion for this particular mode. But interchange part of the mixed 
perturbation does trigger the convection everywhere in the Uy) plane except at Uy = 0 , 
since another condition 7 < 1 -|- a is required for truly mixed mode. It is the convection that 
makes the ‘ridge’ incline with increasing Ux- In the case shown in Figure 9b, the magnetic 
Rayleigh-Taylor instability is activated by the pure interchange mode, as the disk satisfies 
now the condition 7 < 1 — a. This is why the tunnel wall of Figure 9b is partially ‘open’ on 
the ny = 0 side. 

The three-dimensional plots shown in Figures 9c and 9d are for the odd-parity solutions 
but with the same set of system parameters as in Figures 9a and 9b, respectively. We may 
relate the two peaks located on the Ux- and Vy- axes to the Jeans and the Parker-Jeans 
instabilities, respectively. The peak with z/j, = 0 is significantly higher than the one with 
Vy = 0, since the undular mode perturbation provides the Parker instability with favorable 
held configurations to assist the Jeans. On the other hand, the pure interchange mode may 
not get any boosting from the Parker instability. Consequently, the pure Jeans peak on the 
z/a,-axis is always lower than the Parker-Jeans on the Uy-axis. The long ridge along z/y ~ 1 is 
of course owing to the convection-assisted Parker instability; as for the case of even-parity 
solutions, the convection makes the ridge height incline with z/^,. But the ridge reaches a 
hnite level in the limit Ux^oo. 

When mixed mode perturbations of the odd-parity are given to the self-gravitating, 
magnetized, gaseous disk, the Jeans, Parker and convective instabilities may all develop 
in the disk, and their relative importance depends on the (z/^,, z/y) combination. We have 
constructed, from Figure 9c, a plot of equal growth-rate contours in the {i'x,Vy) plane and 
presented the resulting contour map in Figure 10. The thick dashed-line is for a special 
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contour, whose growth rate is equal to the the highest ridge level: Any perturbations with 
wave numbers falling within the dashed-line boundary would grow faster than the ones 
outside. Therefore, with such perturbations, a collaboration of the Parker and the Jeans 
instabilities may well suppress the convection and the disk will bear out large-scale structures. 

In a recent two-dimensional MHD simulation, Kim & Ostriker (2002) followed nonlinear 
responses of a self-gravitating, magnetized, differentially rotating, thin disk to an externally 
imposed spiral potential, and demonstrated regularly spaced spurs jotting out initially in an 
almost perpendicular direction to the arm. Although nonlinear interaction of the Parker- 
Jeans instability and the spiral shock wave is out of the scope of this paper, it should be 
pointed out that the highest point of the growth rate surface is located on the z/^-axis (see 
Figure 10). This is always so, if the disk is thicker than ~0.549i7 (see Appendix C). The 
resulting structure is, therefore, of a cylindrical shape with its long axis being aligned per¬ 
pendicularly to the direction of magnetic field lines or spiral arm (Nagai et ah 1998; Lee 
& Hong 1999; Chou et al. 2000). It is then suggestive that the cylindrical structures of 
this study will make the spurs in the context of differential rotation and density wave. The 
odd-parity undular mode seems to have tailored the regular spacing and initial alignment of 
the spurs. 


4. FORMATION OF LARGE ISM STRUCTURES 

A replacement of driving agent of the Parker instability from the external to the self 
gravity opens up a new route, through which the instability may work constructively with 
the Jeans gravitational instability to form large-scale structures in the Galactic ISM disk. 
In a qualitative sense, the undular mode odd-parity solutions seem to have all the means of 
alleviating the difficulties involved in the time scale, length scale and density enhancement of 
the classical Parker instability. In this section we will extend our qualitative understandings 
of the dispersion relation to making quantitative estimates of the various scales. 

In this study the Galactic ISM disk is characterized by midplane mean density Uq of 
hydrogen nuclei, rms dispersion Cg of “cloudlet” velocities, ratio a of magnetic to gas pressure, 
and thickness ;2a (= CaH) of the disk. Unless tanh(Ca is substantially less than unity, details 
of the dynamics involved with the Parker and the Jeans instabilities do not depend on the 
thickness, which leaves no, Cg and a as the system parameters of principal importance. We 
will examine how these three parameters control the minimum growth time of the Parker- 
Jeans instability, its corresponding length scale, and the mass associated with the length. 
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4.1. A Comparison Standard of the ISM Parameters 

To calculate dynamically important scales of the Parker-Jeans instability, we need to 
establish a ‘standard’ set for the system parameters. By taking averages of various observa¬ 
tions, Ferriere (1995, 1998a, 1998b) presented space-averages of the ISM parameters, which 
may serve us a starting point for setting up the standard. According to his 1998 results, 
the midplane equilibrium number density of hydrogens in the solar neighborhood is Uo = 
1.17 cm“^, which includes all forms of hydrogens residing in the four ISM phases and the 
molecular medium as well. With an inclusion of helium (9% of hydrogen by number), the 
mean molecular weight of ISM per hydrogen nucleus becomes /iisM = 1-36, which makes 
^^ohiSM = 1.59 cm“^. As a comparison standard for density we will simply take Uo/iisM — 
2cm“^ (cf., Holmberg & Flynn 2000; Kulkarni & Heiles 1987). 

Excluding the molecular medium, Ferriere (1998a) determined the total thermal pressure 
of the four ISM phase materials to be 4.81 x 10“^^ ergs cm“^. As the one-dimensional turbu¬ 
lent velocity he took 4.5kms“^ and 6.9kms“^ for the molecular and the cold neutral media, 
respectively, and estimated the total interstellar turbulent pressure to be 6.31 x 10“^^ ergs cm“" 
The sum of thermal and turbulent pressures becomes then 11.12x10“^^ ergscm“^. If his es¬ 
timate of the cosmic ray pressure 9.6x10“^^ ergcm“^ is included in the pressure count, the 
total amounts to 20.7x10“^^ ergs cm“^. With the standard density Uo/iisM — 2cm“^, this 
formally requires the isothermal sound speed or the velocity dispersion to be 7.9 kms“^. This 
seems to be an over-estimate of reality, particularly for ISM near the central plane. Owing to 
the neglect of cosmic rays in this study, we will include only the thermal and the turbulent 
pressures in the isothermal equation of state and take — 5.8kms“^ as a better choice for 
the comparison standard. 

For the magnetic pressure Ferriere (1998a, 1998b) gave 10.3x 10“^^ergscm“^. In terms 
of magnetic field strength, this is equivalent to 5.1/iG, which is signihcantly stronger than 
frequently quoted 3 /xG (cf., Spitzer 1978). If we accept all the pressure values of his estimate, 
a becomes 0.50; if the cosmic ray pressure is not counted, the ratio becomes 0.93. On the 
other hand, if we take 3 /iG for the held strength and count only the thermal and the 
turbulent pressures, then the magnetic-to-gas pressure ratio becomes 0.3. Because many 
uncertainties are involved in deducing the cosmic ray and the magnetic pressures from the 
observed distribution of synchrotron emissivity (Ferriere 1995, 1998a) and because a proper 
treatment of cosmic ray behaviors is beyond the scope of current study (cf., Ryu et ah 2003), 
a may not be hxed at a specihc value at this stage. Instead we will examine how the scales 
vary with a, and chose (uo/iisM — 2cm“^, Cg — 5.8kms“^, 0 < a < 1) as the comparison 
standard. 
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4.2. Time, Length and Mass Scales 


A simple comparison of the odd-parity undular mode solutions to the even-parity ones 
will tell us how much change the constructive interplay between Parker and Jeans instabilities 
can bring to the growth time and the length scale. We have listed in Table 1 the normalized 
maximum growth rates their corresponding wave numbers 

for a selected set of a values. Superscripts “ODD” and “EVN” can be identihed with the 
Parker-Jeans and the Parker instabilities, respectively. The table was prepared with 7 = 1 . 


The ratio given in the fourth column demonstrates that a sizeable reduction is made in 
the growth time by the self gravity. In weakly magnetized disks the time scale is reduced 
to one third the case of pure Parker instability. Since long wavelength perturbations are 
prone to trigger the gravitational instability, we expect the length scale of the Parker-Jeans 
to be longer than that of the Parker instability. Indeed the ratio, Ay,^x) ^he 

last column shows that an about 30 % increase in the length scale can be achieved by the 
Parker-Jeans instability. The relative increase in the length scale varies little with a ; while 
the enhancement in the growth rate changes significantly with the magnetic held strength. 


The minimum growth time, Tgrw = (JJ/cs)|Omax 


tgrw — 3.8 X 10 ’^ years 


1 -|- Q; 

^^ohiSM cm3 


, is given by 
- 1/2 ^ 

|Almax(Qi)| 


(37) 


Because the disk scale height is directly proportional to Cs, the time scale doesn’t depend on 
it explicitly. Substituting the a-dependent D^max Table 1 for the maximum growth rate, 
we examined how the growth time would change with a. Results are shown by thick solid 
line of Figure 11. The two lower curves in thin solid lines are to be identihed with specihc 
cases of rZo/iisM = 2cm“3 and 3cm“3. Please note that the system parameter no/^isM is 
absorbed in the ordinate for the thick solid line. To read the growth time scales for the 
latter two cases directly oh the thin lines, the Uo /iisM factor should be ignored. 


In terms of the system parameters the scale height H becomes 


~ 1.9 X 10^ pc 


1 “h Q; 

1/2 


nofJ-isM cm3_ 


_5 kms~^_ 


(38) 


which serves as the fundamental length for our disk system. The perturbation wavelength, 
A = 27riJz/“3^, corresponding to the maximum growth rate is given by 


A ~ 1.2 kpc 


1 CK 

1/2 

Cs 

1 

.^o/ilSM cm3_ 


_5 kms~^_ 

max(ci) 


(39) 


The observed mean separation of the large-scale structures can be compared directly with 
A, if z/y max is associated with the odd-parity solution. In the case of even-parity, however. 
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the observations should be compared with A/2. When Table 1 is substituted for 

max(tt) in the above equation, the length scale varies with a as in Figure 12, where Uo/iisM 
and Cg are absorbed in the calculation of ordinates for the thick solid line. To read the length 
scales for the two lower curves please ignore the Uq /iisM factor. 

For a mass scale we count total amount of gas material that the initial equilibrium disk 
had within a rectangular column of cross-section A^^max x A^j^max extending vertically from 
—z^ to +Za_- The wave number ^ax of the maximum growth rate hxes azimuthal extent 
of the column cross-section. The undulation of held lines under self gravity introduces to 
the system a natural length scale in the azimuthal direction. However, it doesn’t hx the 
radial extent, since the highest growth rate point in the dispersion-curve surface is always 
located on the z/^-axis. (see, for example. Figure 10.) Only scale we may call in for the 
radial extent is then the critical wavelength of marginal stability or the width of spiral arm. 
In our normalization scheme the critical wave number z/crt becomes 1 or very close to it 
(see Appendix C), which gives a radial extent of Aa;,max — 27rff. Since the total column 
density of halo-bound isothermal gas disk is 2Ffpo(0) tanh(/a, we may have the mass scale, 
Mcrt - 2ifno/iisM Jtih tanh Ca ■ 27rH'z/“^ax ' 27rif, as 


Mcrt — 1-4 X 10^ Mq 


■ Cg 

3 

■ (! + «)=' ■ 

-5kms“^- 


Ao Asm cm3_ 


1 

^y, max (ci 


(40) 


When in Table 1 is substituted for ^ax («) in the above equation, the mass scale 

varies with a as in Figure 13a, where the UoHism factor is absorbed again in the ordinate. 
Since this mass scale depends on the velocity dispersion most sensitively, we have illustrated, 
in the figure, for the following four cases specihcally: Cg = 1, 2, 3, and 5 kms“^. 


Since the disk scale height is about 200 pc, the critical wavelength often turns out larger 
than the width of spiral arm, which is probably about a half kpc. For such cases, the arm 
width, AR, should be taken for the radial extent of the rectangular column. The resulting 
mass scale. Alarm, now reads 


Alarm ^ 1-2 X 10^ Mq 


Cg 

2 

'AR 


1 -|- o; 

-5kms“^- 


kpc 


_^y, max (ci) _ 


(41) 


It is interesting to notice that once the arm width is hxed by the global dynamics of rotating 
disk. Alarm docsu’t depend on the midplane density explicitly. Substituting for the 

wave number of the maximum growth rate, we have calculated the mass scale as a function of 
the magnetic-to-gas pressure ratio, and shown the results in Figure 13b for the four selected 
cases of the velocity dispersion. This hgure was prepared with AR = 1 kpc. 
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4.3. GMCs and HI Superclouds 


The two most prominent types of large interstellar structures are the GMCs (Dame 
et ah 1986; Blitz 1993) and the HI superclouds (Elmegreen & Elmegreen 1987; Efremov 
& Sitinik 1988; Elmegreen 1994). Clouds of each type show substantial spreads in mass, 
size, mean separation, internal velocity dispersion, and density. Particularly in the velocity 
dispersion the two types almost overlap with each other (Falgarone & Lequeux 1973; Larson 
1981; Heiles & Troland 2003). However, the GMCs differentiate themselves as a distinct type 
from the HI superclouds in most of the directly observable properties (Alfaro et al. 1992). 
In Table 2 we have summarized observations of GMCs and HI superclouds separately. For 
entries of the GMCs we relied heavily on Blitz (1993), and the HI supercloud properties 
are mostly from Elmegreen & Elmegreen (1987). The mean separation of the superclouds 
is known to depend on mass: The heavier are the clouds, the wider becomes their mean 
separation. 


It is not easy to have a relevant time scale for structure formation directly from obser¬ 
vations, but a number of estimates can still be made. For example, time required to pass an 
arm will serve a good criterion for the growth time, since gathering activity for the structure 
formation should take place mostly while the gas material stays in a spiral arm. As suggested 
by Kim & Ostriker (2002), one half of the arm-to-arm moving time can be a reasonable mea¬ 
sure for the arm-crossing period. In an m-armed spiral galaxy it will take 27ri?/[mi?(r2 —Dp)] 
for the gas material to complete a journey along the galactocentric circle of radius R from 
one arm to the next. Here, D is the angular velocity of the galactic rotation at galactocentric 
distance R and Dp is that of the spiral pattern. With an approximation Dp ~ D/2 (Binney 
& Tremaine 1987) the arm-passing time simply becomes 47r/(mD). Keeping in mind the 
solar neighborhood in the Galaxy, we take 220kms“^ for the rotation velocity at i?=8.5kpc, 
and to the number of arms we may assign four (cf. Vallee 1995) to have 


tapt ^ 60 X 10® years 


(-) 

26kms ^ kpc ^ 

\m J 

D 


(42) 


Larson (1994) estimated the CMC growth time from the gas consumption rate in the 
Galaxy. The star formation rate is in the order of 3 Mq per year, and about 2 % of cloud 
mass is known to be converted into stars (Myers et al. 1986; Evans & Lada 1991). This 
means that 150 Mq of gas is being condensed, in every year, into star-forming molecular 
clouds. Since the total amount of gas in the Galaxy is about hxlO^M©, mostly in the form 
of GMGs, the time required to collect this much gas into the GMGs is about 33x10® years. 
One may then take 30x10® years as the formation time scale of the GMGs. 

We now have three time-scale values: The arm-to-arm crossing period, 1.2x10® years. 
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is thought be a somewhat generous time-scale criterion for the formation of HI superclouds, 
while the arm-passing time, 60x10® years, is a rather strict one for them. On the other hand, 
the gas consumption period, 30x10® years, will serve as a reasonable time-scale criterion for 
the formation of GMCs. 

As is demonstrated in Figures 11, 12 and 13, the ranges of time, length and mass scales 
that are expected from the undular odd-parity solutions compare well with the observational 
criteria of the HI superclouds summarized in Table 2. To be specihc, with the system param¬ 
eter set, Uo/iisM — 2cm“^, Cg ~ 5.8kms“^ and a ~ 0.5, we have the growth time tgrw — 
4.3x10'^years, the mean separation A ~ 1.9kpc, and the mass scale Alarm — 2.8xlO^M0. 
With the same set, the disk scale height becomes H ~ 190 pc, and the resulting critical 
wavelength, 1.2 kpc, extends in the radial direction beyond the spiral arm. As illustrated in 
the hgures, the choice of a = 0.5 won’t make a critical issue in bounding the various scales; 
the scales depend least sensitively on a, unless it is much bigger than unity. From the com¬ 
parison we conclude that under canonical conditions of the Galactic ISM the Parker-Jeans 
instability can drive the ISM disk to form condensations of the HI supercloud scale. 

To make the theoretical scales agree with the observational criteria of the GMGs, one 
has to strain the system parameters to extreme degrees. For example, a trial set of Uq /xism — 
50cm“^, Cs — 1.2kms“^ and a ~ 1.0 gives us 11 million years for the growth time, 510 pc for 
the mean separation, and Alert — 1.4xlO®M0 for the mass scale of the GMGs. As done in 
this example, one could devise a combination of the system parameters that would formally 
satisfy the observational criteria of the GMGs. But the choice Uq /iisM — 50 cm“^ seems 
too high to be the density of initial equilibrium conhguration. And the choice of this high 
value of density requires the velocity dispersion of cloudlets too low. The undular odd-parity 
solutions under the self gravity alone may not grow the GMG and the HI supercloud scale 
structures simultaneously. 


5. DISCUSSION AND CONCLUSION 

Through the present study we come to realize that in previous investigations role of the 
interchange mode was unduly stressed in its tendency of disrupting large-scale ISM structures 
that the undular mode might otherwise have formed. The unfair emphasis was the price to 
pay for both adoption of the external gravity as and neglect of the self gravity from driving 
agent of the Parker instability. In this paper we have examined what modiheations could be 
made to the classical picture of the Parker instability by assigning the driving task solely to 
the self gravity. 



If the vertical acceleration rendered by the adopted driving agent is a constant of per¬ 
pendicular distance from the central plane, it may act as a source of buoyancy force over 
the entire disk. Consequently, under such unrealistic uniform external gravity, convective 
motions are to be activated everywhere in the disk, and the interchange mode will indeed 
dominate the magnetized gas disk over the undular one. In reality, however, the vertical 
acceleration due mostly to the Galactic stars increases linearly, from zero at the midplane, 
with the perpendicular distance up to about a half kpc, beyond which it remains more or less 
constant. This type of external gravity can not be a buoyancy source in the disk midplane at 
least, where GMCs are mostly found; while the convection is pronouncedly active, under the 
realistic non-uniform gravity, only in high altitude regions, where the vertical acceleration 
reaches its maximum. 

The convection wouldn’t then be an obstacle for the classical Parker instability to form, 
for example, GMGs in the midplane, if the undular mode could collect enough material 
there before the interchange mode develops chaotic sheets at high altitudes. However, as 
long as the buoyancy force is available, the growth rate of the convection increases within the 
limiting value as the perturbation wavelength gets smaller and smaller. Therefore, even if all 
the realistic features are taken into account for the external gravity, the convection triggered 
by the interchange mode would develop small-scale features well before the classical Parker 
instability excited by the undular mode generates any dense large-scale structures in the disk 
central plane. Magnetic valleys eventually anchor some material towards the central plane, 
but the material takes form of vertically oriented thin sheets rather than of dense clumps, 
because the external gravity has only vertical component of acceleration. As long as the 
driving task is assigned to the external gravity, there is no means to divert the convective 
motions from the vertical direction and to form clumps instead of sheets. 

The thick dashed-line in Figure 10 represents the growth rate contour corresponding to 
the highest ridge level. Any perturbations whose wave numbers are within the dashed-line 
boundary would grow faster than the ones outside. Therefore, with such perturbations, a 
collaboration of the Parker and Jeans instabilities may well suppress the convection and the 
disk will bear out large-scale structures. 

Once the self gravity takes over the driving task from the external gravity, the effect of 
the Jeans gravitational instability due to the self gravity appears on the dispersion curves. 
The maximum growth rate of the Jeans instability at the large-wavelength perturbations 
can be be higher than the limiting growth rate of the convection at the small-wavelength 
perturbations, depending on the adiabatic index (cf.. Figure 7). This is not what we’ve seen 
from the external gravity case. Therefore, for long-wavelength perturbations of odd-parity, 
the maximum growth rate of the undular mode becomes even faster than the highest growth 



rate value of the interchange mode (cf., Figure 10). The odd-parity undular mode thus 
opens up a route in the magnetized disk, through which the Parker and Jeans instabilities 
work jointly to form large-scale condensations in the midplane before the interchange mode 
disperses them into thin sheets. 

In this way the Parker-Jeans instability is expected to remove major difficulties in form¬ 
ing the large-scale ISM structures. To check whether this expectation can be fulhlled by the 
odd-parity undular mode, we have compared the resulting time, length and mass scales with 
the corresponding observations of the HI superclouds and GMCs. 

As far as the time scale is concerned, the Parker-Jeans instability grows quickly enough 
to satisfy the observational criteria for the HI superclouds and the GMCs as well (cf.. Figure 
11). The formation time scales inferred from various observations don’t impose any tight 
limits on the system parameters, particularly on a, and the comparison done in Figure 11 
doesn’t seem to differentiate the HI superclouds from the GMCs either. It was a bit surprise 
to us. Because of the neglect of cosmic rays, we thought, the time scale of the Parker-Jeans 
would meet the observational criteria only barely. In a recent study Ryu et al. (2003) treated 
cosmic rays with the diffusion-convection equation and realistic diffusion coefficients for both 
directions, perpendicular and parallel to the held lines. They have demonstrated that by 
pnmping an extra buoyancy into the disk system, cosmic rays make an abont factor two 
increase in the growth rate of the classical pure Parker mode, i.e., of the instability triggered 
by the nndnlar pertnrbation nnder the nniform external gravity. The factor two may not 
be fully applied to the Parker-Jeans case, because the uniform external gravity exaggerates 
bnoyancy. Nevertheless, if cosmic rays had been included, a substantial farther redaction 
would have been made in the growth time scale. It may worth reminding that the Jeans 
gravitational instability was called in to snpport the Parker in the hrst place. And the resnlt 
was snccessfnl, as was demonstrated by the |f^y,^xl/l^y,^xl ” ratios in Table 1: A small 
increase in will result in a sizeable amplihcation throngh the ratio. An inclusion of 

rotation wouldn’t change the sitnation in any signihcant way, because the rotation is known 
to decrease the growth rate by a few to hfteen percents only (Shn 1974). Along this line of 
reasonings, if 2 cm“^ is not an nnrealistic choice for Uofiisu in the Galaxy, we believe that 
the time scale we now have is an under-estimate for the Galactic ISM disk. The reason for 
the nnderestimation seems to be in onr neglect of the external gravity. 

As to the length scale, the Parker-Jeans instability seems to operate more efficiently 
than what is required for the HI superclouds. Under ‘canonical’ conditions of the Galac¬ 
tic ISM with a = 0.5, the odd-parity nndnlar mode shows its maximum growth rate at 
wavelength 1.5 kpc, which is less than the observed mean separation, 2.4 kpc, of the HI 
superclonds (McGee & Milton 1964; Elmegreen & Elmegreen 1983; Alfaro et al. 1992; Efre- 
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mov 1995). Because of the conceptual difficulty in applying a single velocity dispersion to 
the entire ISM disk, one may not consider this much difference in the length scale to be 
serious. However, we think the difference is a signihcant one, because existence of cosmic 
rays would reduce the wavelength further. Therefore, comparison of the length scale also 
suggests an over-subscription of the self gravity in our disk system. 

From dispersion relations shown in Figures 1 and 2 we notice the growth rate peak 
for the even-parity undular mode doesn’t survive in the competition with the odd-parity: 
The growth rate peak of the pure Parker mode is overwhelmed by the Parker-Jeans peak. 
This holds true for wide ranges of disk thickness and magnetic-to-gas pressure ratio. As 
was illustrated in Figure 3, the Parker peak may survive the Parker-Jeans overwhelming 
only in such unrealistically soft media as the ones having 7 <0.7, for example. In the 
magnetized gas disks where only the self gravity is operative to excite instabilities with the 
undular mode, it doesn’t seem to be possible to have two distinct growth-rate peaks in 
the dispersion curve, and the present study indicates that the overwhelming Parker-Jeans 
maximum is strongly tied to the structures of the HI supercloud scale, not to the GMC 
scale ones. As Elmegreen(1994) pointed out, the characteristic mass and mean separation 
of the HI superclouds are natural outcomes of the gravitational instability in the self gravity 
dominating disk systems. 

To relate both the GMC and HI supercloud structures to the single Parker-Jeans peak, 
one may rely on two-stage formation scenarios: The HI supercloud forms hrst and some time 
later the GMCs emerge as fragment products of the supercloud. This idea was suggested in 
the previous section reluctantly, because of the GMCs that are residing outside the super¬ 
cloud. Almost all of the HI superclouds contain GMCs, but not all the GMCs belong to the 
HI superclouds (Elmegreen & Elmegreen 1987; Efremov & Sitnik 1988). Instead rehning the 
two-stage idea, we should hrst pay attention to the external gravity heretofore ignored. In 
fact the external gravity renders vertical acceleration about 5 times stronger than the one ex¬ 
pected from the equilibrium ISM stratihcation in our Galaxy (Bienayme et ah 1987; Ferriere 
1998a; Boulares & Cox 1990; Lee et al. 2004). Of course in later stages the ever increasing 
self gravity of ISM eventually override the hxed external gravity, but in the initial stage the 
self gravity is dehnitely over-subscribed in our model disk. A balanced subscription between 
the external and self gravities for the task of instability driving would reduce dominance of 
the Parker-Jeans peak on one hand and enhance the Parker maximum on the other. In this 
way the odd-parity undular mode may exhibit two distinctly separated growth-rate peaks in 
the dispersion curve, with the one at smaller wave number corresponding to HI superclouds 
and the other at larger number to GMCs. 


From these discussions we may conclude the paper as follows: The present study has 
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reconfirmed that the Parker instability triggered by the undular mode under external grav¬ 
ity may not be a route towards the formation of such large-scale structures as the GMCs in 
the Galaxy (Kim et ah 2004). Even under the self gravity, the even-parity undular mode 
solution can not accommodate important observational properties of the large-scale struc¬ 
tures; however, the odd-parity undular mode is shown to have all the potentials to remove 
major difficulties that have been confronting the classical Parker instability as the formation 
mechanism of the large-scale structures. The resulting time, length, and mass scales of the 
Parker-Jeans instability are shown to be in better agreement with the HI superclouds than 
with the GMGs. An inclusion of the external gravity along with a proper treatment of the 
cosmic rays would rehne the agreement towards the desired direction. 


It is our pleasure to acknowledge illuminating discussions and stimulating conversations 
with W.-T. Kim, B. G. Elmegreen, J. Kim, D. Ryu, and J. Franco. The authors are grateful 
to B.-G. Koo and T. Nakagawa for guiding us to some important literatures. We are also 
thankful to J. Pyo for his assistance in the preparation of some of the diagrams. While prepar¬ 
ing this paper, SSH was benehted greatly from Astronomical Information Research Group 
of the KASI, Korea and also from Laboratory Infrared Astrophysics of the ISAS/JAXA, 
Japan. This study was supported in part by KOSEF through grant R14-2002-058-01003-0. 



27 


APPENDIX 


A. Coefficients of Perturbation Equations 


The coeffients of the two second order ODEs [Eqs. (24) and (25)] are as follow: 
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Here k‘^ means and the other symbols have the same meanings as in §2 of the main 

text. 


B. Stable Modes 

In self-gravitating, magnetized, gaseous disks, we expect to have such stable wave modes 
as gravity wave, Alfven wave, and fast and slow MHD waves. To check whether our system 
of the linearized perturbation equations indeed has all these modes, we take the equations 
in various limits of perturbation wave numbers. 



B.l. Vertical Gravity Wave Mode 


The vertical gravity mode can be isolated by taking the perturbation equations in 
the limit i^x = ^ and Uy = 0. This wave travels vertically and is related to convective 
motions (Lamb 1945). 

We eliminate ip from equations (24) and (25) and put = 0 and = 0 in the resulting 
equation for Uz to have 
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hold true, we may integrate equation (Bl) and have 


^ d 1 -|- a 

-j^Uz - ^Q^Uz - -—^^ Uz = 0 . 

dC dC 2a + 'y 

Changing variable ( to Z = tanh(C, we rewrites this as 
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where a is dehned as 
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We then substitute H = {1 — Z‘^)~Puz{Z) ioi Uz{Z) andw = (l-|-Z)/2 for Z. With parameter 
p dehned by 

= 4p(p-1/2), (B6) 

Equation (B4) reads 

W2 7 

vj{l — vj) — 2 -^ + + 1/2 — (4p -|- l)vj} -— D — Ap^D = 0 , (B7) 

dw dzu 

which is a hypergeometric differential equation. In terms of Z its general solution is given 
by 

Uz{Z) = {1 - Zy [C,F {2p, 2p-, 2P+1/2-,1/2(1+ Z)} 

+ Cs {1/2(1+ Z)}B2-2Pi7’{1/2,1/2; 3/2-2p; 1/2(1 + Z)} 1 , (B8) 
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where F is a hypergeometric function, and Ci and C 2 are integral constants (Polyanin & 
Zaitsev 1995). 

From the boundary conditions (Eq. [B2]), we construct a set of two simultaneous indef¬ 
inite equations of Ci and 6 * 2 • To have nontrivial solutions for the integral constants, matrix 
of the simultaneous equations should have vanishing determinant G{p). This condition gives 
us secular equation G{p) = 0 with the following dehnition of G{p): 

tP' ( ‘2t) “t“ li ^ 

= {2p + l/2){2p + 3/2) ~ 

r o 2 p-l /2 

- (V - 1/4) F(1/2, 1/2; 3/2 - 2p; 1/2) + F(3/2, 3/2; 5/2 - 2p; 1/2) (B9) 

+ 8(4p - 3)(4p - 5) ^^^>1 ''PP-■ 

In Figure A 1 we have plotted G{p) against p. As can be seen from the fignre, G{p) 
becomes zero at negative half integers and also at zero. In terms of p the resulting eigen- 
frequencies fl are given as 

n^ = ^^^2p{2p-l), (BIO) 

1 -|- q; 

where p= - 3k{2k- l)/2 with A; = 0, 1, 2, ■ ■ 

This mode of stable waves represents gravity waves propagating vertically in the eqni- 
librinm disk ( cf., Eq. [50] in Shn 1974). In the absence of magnetic helds Lamb (1945) called 
this as the atmospheric wave. This mode appears at the very beginning stage of nonlinear 
MHD simulations. 


B.2. Fast MHD Wave 

For the limiting case ^ p 0 and p = 0, perturbations in the x—direction disturb density, 
pressnre and magnetic fields in such a way that the field strength and density show alternating 
regions of compression and rarefaction. This is a signatnre of simple acoustic waves. However, 
it is not easy to isolate this type of acoustic waves near the dense midplane, where the self- 
gravitational force is likely to override the restoring force of the magnetic pressnre. In order 
to isolate the mode clearly, we apply the eqnation to the limit z^oo, where the gravity force 
becomes very weak. 

In this limit p 0, p = 0, z ^ 00 ) one may easily isolate the fast MHD wave mode, 

1 + a ^ ’ 
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or 

u^ = {vl + al)kl (B12) 

where al means 70^. The numerator 2 q ; + 7 of equation (B12) indicates a coupling of Alfven 
and sound waves. And the ,^^-term suggests that the wave propagates in the horizontal 
direction perpendicularly to the initial magnetic helds. The denominator 1 + a comes from 
the dehnition of the scale height, H. 


B.3. Alfven and Slow MHD Waves 


To isolate the Alfven and the slow MHD modes, it is better to suppress the Jeans and 
the Parker instabilities. We thus push the perturbation equations to the limit 1 7 

and have the following dispersion relation: 
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and the slow mode of the MHD waves 
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The two waves propagate in the horizontal plane along the direction of initial magnetic helds. 
The 20 : 7 /( 20 ; + 7 )-term represents the cusp speed for the group velocity of the slow wave. 


Finally, the magnetized sound wave is derived from the y-component of the momentum 
equation. The sound waves propagate along the magnetic held direction, with the dispersion 
relation given by 




—— 7 ^ , or = Qs^ky^. 
1 + 0 '’ ^ 


(B16) 


C. Marginal Stability Analysis 

For thick disks the Parker instability may assist the Jeans gravitational instability, so 
the maximum of the growth rate curve is located on the z/^-axis in the plane of perturbation 
wave numbers. This was demonstrated in Figures 9c, 9d and 10, where the peak on the 
z/y-axis is higher than the one on the z/^j-axis. As the disk thickness decreases, the growth 
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rate peak on the z/^-axis decreases, since the Parker instability is harder to get triggered in 
thinner disks. In the disk of the critical thickness, the Parker instability disappears, and 
the two peaks on the and z/^^-axes become of the same height. For the disk thinner than 
the critical thickness, the magnetic helds obstruct the system from developing the Jeans 
instability along the y-direction. So, the maximum of the growth rate curve is now located 
on the z/a;-axis. 


As was shown in Figure 3, the dispersion relation of the Jeans-Parker instability has 
two local maxima. The peak at smaller wave number is for the Jeans instability, while the 
one at larger number for the Parker instability. The dispersion curve crosses a critical wave 
number, which separates the region of stability from that of instability. The critical wave 
number comes into being due to magnetic tension, and the Parker instability gets stabilized 
by the tension. Thus, investigating the marginal state, we can hnd the critical thickness for 
the disk that may trigger the Parker instability. 


With .^ = 0 and hi = 0, equation (21) describes the marginal state of the Parker instability 


as 
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where all the symbols have their usual meanings. Applying equation (Cl) to the disk bound¬ 
ary and denoting the eigenvalue of the equation by 
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where ©a = tanh , we have the critical wavelength as 
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Because the critical wavelength, Acrt.non, of the non-magnetized, isothermal, gaseous disk 
which is inhnitely extended in vertical direction is given by 
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we know, equation (C3) should reduce to Equation (C4), in the limit a —> 0, 7 = 1 and ©« 
= 1. This can be realized by simply assigning unity to EV. From this we infer EV = 1 is 
an eigenvalue of our halo-bound hnite disk. Therefore, we may have the normalized critical 
wave number of the perturbation propagating along the magnetic helds as 

[ 2(1g - 7)(1-F g) 2 

Vert = - ©a - 1 • (C5) 


gy 



32 


This nicely recovers the critical wave number, \^2a + 1, of a magnetized, isothermal, gaseous 
disk infinitely extended in vertical direction (Stodolkiewicz 1963). 


Since the Parker instability may not be triggered in a disk whose critical wave number 
is zero, from equation (C5), we have the critical adiabatic index for an onset of the Parker 
instability as: 


bcrt — 


02 + T(i + 202)a 


(C 6 ) 


which again recovers the well-known = (1 + ck:)^/(l + 3a/2) condition for the classical 
Parker instability in an infinitely extended disk under external gravity (Kim & Hong 1998). 
In the disk with 7 < 'jcrt, we won’t expect the Parker instability to take a place. In Figure 3 
we have illustrated how the material rigidity controls the dispersion relation of the Parker- 
Jeans instability. The dispersion curve usually shows two bumps, one from the Jeans and 
the other the Parker. But the Parker bump becomes inconspicuous as the adiabatic index 
increase, the even-parity dispersion curve in dashed line almost coincides with the wave 
number axis when 7 = 1.04, and completely disappears as 7 further increases; from equation 
(C 6 ) we indeed have ycrt = 1-05 with a = 0.1 and Ca = 5. 


In §3.1, we have noticed that the Parker instability may assist the Jeans instability only 
in thick enough disks; thin disks may not conceive the Parker instability in the first place. 
What is then the critical thickness separating the thin from thick disk? This can be answered 
from equation (C5) by finding out such 0^ (=tanh(Ca) that may make ^crt zero, namely, 

11/2 

= tanh“^ I -r . (C7) 
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When C^a < Ca.crt, the Parker instability cannot operate in the system, and the magnetic 
fields instead obstruct the self-gravitational instability. Maximum of the growth rate curve 
for the gravitational instability is located, in this case, on the u^-axis, so the gravitational 
instability would generate a long cylinder whose axis is parallel to the magnetic field lines. 
On the other hand, when (a > Ca,crt, the Parker instability assists the gravitational instability, 
and the Parker-Jeans instability is operative in the system. Maximum of the growth rate 
curve for the Jeans-Parker instability is placed now on the z/^-axis, so the instability would 
align the cylinder axis perpendicularly to the field lines. When (a = Ca,crt, the instability 
takes place, and the Parker instability cannot assist the Jeans to grow faster along the 
fields. Thereby, the perturbation propagating along the y-direction coordinates with that 
propagating along the x-direction. The critical thickness, when a = 1.0 and 7 = 1.0 for 
example, becomes 0.549 H. As shown in Figure A 2 , the undular mode {fx = 0) dispersion 
curve in solid line looks almost the same as the interchange mode (z/^ = 0 ) one in dashed 
line. 
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Table 1. Comparison of the Odd- and Even-Parity Solutions 


a 


IqEVN I 

1 ^y,max 1 

ODD /^EVN 
' min / ' min 

,.ODD 

^y,max 

EVN 

y,max 

\ODD /\EVN 
^y,max/ ^y,max 

0.00 

0.665 

_ 

_ 

0.476 

_ 

_ 

0.25 

0.722 

0.259 

0.359 

0.547 

0.750 

1.371 

0.50 

0.766 

0.372 

0.486 

0.621 

0.838 

1.350 

1.00 

0.835 

0.513 

0.614 

0.749 

1.000 

1.334 

2.00 

0.924 

0.675 

0.730 

0.943 

1.239 

1.313 

3.00 

0.981 

0.765 

0.780 

1.093 

1.417 

1.296 





Table 2. 


Observational Properties of GMCs and HI superclouds 


Properties 

GMCs 

HI Superclouds 

Mass [Mq] 

(l~10)xl0'''’ 

(l~40)xl0®'= 

Mean Separation [kpc] 

0.4 ~ 0.6'' 

~ 2.4® 

Velocity Dispersion [km s“'] 

0.4 ~ 8.5'’g 

1.1 ~ 8.P 

Mean Density [cm“^] 

~ 10^ Ha'' 

>9H® 

Size [pc] 

~ 45'' 

~ 400® 

formation time scale [years] 

30x10® 

80x10® 


^arson 1981 

'’Dame et al. 1986 

'’Elmegreen & Elmegreen 1987 

''Blitz 1993 

®Alfaro et al. 1992 

'Heiles & Troland 2003 

^Falgarone & Lequeux 1973 
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Fig. 1.— Dispersion relation of the undular mode = 0) for (a) non-magnetized and (6) 
magnetized gaseous discs of various thickness. The abscissa denotes effective wave number 
Vy, and the ordinate does square of the growth rate. Some of the disk and perturbation 
parameters are given in the upper left corner of each frame. The solid and dashed lines in 
{h) are for odd- and even-parity solutions, respectively. Because even-parity can not activate 
gravitational instability, there is marked suppression in the maximum growth rate from the 
odd- to the even-parity solutions. 
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Fig. 2.— Dispersion relation of the undular mode {fx = 0) for various values of a. The 
abscissa denotes effective wave number and the ordinate does square of the growth rate. 
Some of the disk and perturbation parameters are given in the upper left corner. The solid 
and dashed lines are for odd- and even-parity solutions, respectively. 
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Fig. 3.— Dispersion relations for selected values of 7 . The abscissa denotes effective wave 
number Vy and the ordinate does square of the normalized growth rate of undular mode 
perturbation {y^ = 0.0). Solid and dashed lines are for the odd- and even-parity solutions, 
respectively. The disk parameters are specihed in the upper left corner, and 7 -values are 
marked to each curve. 
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Fig. 4.— Eigen-solutions of the undular mode with odd-parity (a, b) and with even-parity 
(c). Disk and perturbation parameters are specified on top of each panels on the left hand 
side. The abscissa and ordinate are y and 2 ; coordinates, respectively. The left hand side 
panels illustrate the density held, and color changes from dark red to white as logarithm 
of pi{y,z)/po{z = 0) increases. Arrow vectors in the right hand side panels illustrate the 
velocity held with magnitude reference being given on top of each frame. The thin solid lines 
represent the magnetic held lines and the thick lines do the disk-halo boundary. Case (a) is 
an example of the self-gravity dominating Parker-Jeans instability, while in (6) convection 
dominates, particularly in high altitudes. In case (c) the gravitational instability is not 
triggered, because even-parity condition is imposed upon the solution. 
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Fig. 5.— Dispersion relation of the interchange mode {fy = 0) for various disk thickness. The 
abscissa denotes effective wave number and the ordinate does the growth rate squared. 
To each curve is marked half thickness of the disk, and other disk parameters are given in 
the upper left corner. 
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Fig. 6.— Dispersion relation of the interchange mode {uy = 0) for a number of selected 
a’s. The abscissa denotes effective wave number and the ordinate does the growth rate 
squared. To each curve is marked a value, and other disk parameters are given in the upper 
left corner. 



45 


C\? 

a 



7 

0.68 


0.72 


0.76 


0.80 


0.84 


0.88 


Fig. 7.— Dispersion relation of the interchange mode (z/ = 0) for selected 7 values. The 
abscissa denotes effective wave number and the ordinate does the growth rate squared. 
To each curve 7 value is marked, and other disk parameters are specified in the upper left 
corner. The numbers with arrow on the right hand side of the frame indicate square of the 
maximum growth rate reached in the limit > cx). 
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Fig. 8.— Eigen-solutions of the interchange mode with odd-parity (a) and with even-parity 
{b). Disk and perturbation parameters are specified on top of each panels on the left hand 
side. The abscissa and ordinate are x and 2 ; coordinates, respectively. The left hand side 
panels illustrate the density held, and color changes from dark red to white as logarithm 
of pi{y,z)/po{z = 0) increases. Arrow vectors in the right hand side panels illustrate the 
velocity held with magnitude reference being given on top of each frame. The thick lines 
represent the disk-halo boundary. Case (a) is an example of the Jeans gravitational instability 
triggered by long wavelength perturbation, while in (b) the convective motions activated by 
short wavelength perturbation are apparent, particulary in high altitudes. 
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Fig. 9.— A three-dimensional surface plot for the dispersion relation of the mixed mode. In 
all four cases the halo is placed at (a = 5. Other parameters are as follow: In case (a) a = 
1 and 7 = 1 with even-parity; (b) a = 0.1 and 7= 0.8 with even-parity; (c) a = 1 and 7 = 
1 with odd-parity; and (d) a = 0.1 and 7 = 0.8 with odd-parity. The ordinate represents 
square of the growth rate and the two abscissae denote the horizontal perturbation wave 
numbers. In (a) and (b) the gravitational instability has not been activated; while in (c) 
and (d) the Parker-Jeans instability is operative. On the other hand, in (b) and (d), the 
convection criterion 7 < 1 — a for the pure interchange mode is fulhlled, and hence the 
Rayleigh-Taylor instability gets triggered. The inclining ‘ridge’ of the growth-rate surface 
reaches in the ^ 00 a hnite level, which is lower than the Parker-Jeans peak on the i/y 


axis. 
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Fig. 10.— Equal growth rate contour map. This is constructed from the 3-dimensional 
surface plot of the dispersion relation for undular odd-parity solutions with = 5.0, a = 1.0, 
7 = 1.0. The dashed line contour in blue corresponds to the highest ridge level. Therefore, 
the perturbations whose wave number pair falls within the blue contour are expected to grow 
faster than the convective instability. 
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Fig. 11.— Growth time scale versus magnetic-to-gas pressure ratio. Notice that the Uo/iisM 
factor is absorbed in the ordinate. When the ordinates for the lower two curves in thin solid 
line are read directly off the hgure, the density factor should be ignored. 
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Fig. 12.— Wavelength of maximum growth rate versus magnetic-to-gas pressure ratio. No¬ 
tice that Uo/iisM and velocity dispersion are absorbed in the ordinate for thick solid line. 
When the ordinates for the lower two curves in thin solid line are read, the multiple of the 
two factors should be considered unity. 
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Fig. 13.— Mass scales based on the critical wave number (a) and the spiral arm width (b) are 
plotted against magnetic-to-gas pressure ratio for four selected cases of velocity dispersion. 
Notice that the no/^isM factor is absorbed in the ordinate of (a). Once the global dynamics 
of rotating, self-gravitating, magnetized, gaseous disk hxes the arm width for us, the mass 
scale based on the arm width doesn’t depend on the midplane density explicitly; it does on 
velocity dispersion to the second power. To prepare this figure AR = 1 kpc is used. 




























53 



Fig. A2.— Dispersion relations of the critically thick disk. The abscissa represents the 
normalized wave number Vy for the undular mode (solid line), and for the interchange 
mode (dotted); while the ordinate denotes square of the normalized growth rate. Some of 
the parameters used in the calculation are given in the upper left corner of the frame. 




